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Abstract 


In  turbulent  reacting  flows,  completion  of 
chemical  reactions  requires  diffusion  of  reactive 
species  at  a  molecular  level.  Eddies  rich  in  one 
reactant  can  be  present  in  many  turbulent  flows 
involving  finite  rate  chemical  reactions.  Local 
species  stratifications  are  dissipated  by 
molecular  diffusion  and  chemical  reaction  under 
the  influence  of  the  local  flow  field.  In  many 
turbulent  flows,  the  local  flow  field  has  a  strong 
rotational  character.  It  is  mostly  this 
rotational  motion  which  enhances  mixing  through 
flow  stretch  and  convection  leading  to  mixing  at 
the  molecular  level. 

In  this  paper,  molecular  mixing  and  finite 
rate  chemical  reactions  in  a  two  dimensional 
viscous  vortex  are  examined  analytically.  Two 
species  initially  separated  across  a  plane  are 
allowed  to  diffuse  and  react  in  the  presence  of  a 
line  vortex  centered  at  this  separation  plane. 
Concentration  fields  and  the  probability  density 
functions  (pdf)  are  determined  for  a  range  of 
vortex  strengths  and  for  several  values  of  Schmidt 
number  at  different  times  during  the  growth  of  the 
vortex.  Finite  rate  isothermal  chemical  reactions 
are  also  considered  The  results  for  the  reacting 
vortex  are  compared  with  those  for  the  nonreacting 
vortex  and  some  insight  is  gained  concerning  the 
form  of  the  pdf  in  such  configurations. 


Method  of  Analysis 


Nonreacting  Vortex 

The  concentration  field  within  a  nonreacting 
viscous  line  vortex  Imposed  coincident  with  an 
interface  separating  the  two  species  is  calculated 
by  analyzing  the  local  behavior  In  a  Lagranglan 
frame  of  reference  Computation  of  the  successive 
locations  of  a  material  element  starting  from  its 
Initial  position  is  achieved  by  integration  of  the 
flow  field  with  respect  to  time.  In  this  study, 
the  analytical  solution  for  the  flowfield  of  a 
viscous  line  vortex  is  utilized.  Across  the 
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material  surface  separating  the  two  species, 
molecular  diffusion  takes  place  as  this  interface 
is  stretched  and  transported  in  the  vortex  flow 
field  The  local  concentration  distribution 
around  interfaces  can  be  obtained  analytically. 
Until  this  point,  the  analysis  generally  follows 
References  2  and  3.  In  the  solution  for  the 
concentration  field  within  the  vortex,  the 
concentration  at  a  given  location  is  composed  from 
superposition  of  analytical  solutions. 
Concentrations  at  many  locations  are  obtained  and 
used  to  construct  the  pdfs  assuming  periodicity  of 
the  vortices  and  utilizing  Taylor's  frozen  flow 
assumption . 


Chemically  Reacting  Vortex 

In  case  of  finite  rate  chemical  reactions 
within  the  vortex,  an  approximate  method  based  on 
Green's  function  solution  of  the  reacting  species 
conservation  equation  is  implemented  Solution 
for  the  mass  fraction  of  one  reactive  species, 
together  with  the  solution  for  the  appropriate 
Shvab- Zel ' dovich  variable,  allows  determination  of 
both  fuel  and  oxidizer  distributions  within  the 
vortex . 

In  this  paper,  a  dilute  isothermal  reaction 
between  NO  and  0-^  is  considered.  The  reaction 
products  are  NO 


and  O2  and  both  reactants  are 


diluted  in  Nj .  The  mass  fractions  far  from  the 


vortex  are  taken  to  be  0.01  for  both  reactants. 
The  second  order  rea.  tion  rate  was  taken  from 
Givi,  et  al  ^  ^  . 

Detailed  description  of  the  analyses  for 
nonreacting  and  chemicslly  reacting  vortices  can 
be  found  in  Reference  (5). 


Results 


The  first  four  figures  present  results  for 
the  nonreacting  case.  Of  course,  they  also  apply 
to  the  passive  scalar  in  the  reacting  case 
Figure  la  shows  the  distortion  of  an  initially 
horizontal  material  line  located  along  the  x  axis 
at  a  certain  instant  after  the  vortex  initiation 
for  a  nondimensional  vortex  strength 

(Re  -  r/2*v)  of  50.  The  initial  concentration 
field  corresponds  to  the  upper  half  plane  of 
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cor.,  t* t  a:  ior.  0  0  and  the  lower  half  plane  of 
cor.,  er.rra:  ipr  1  0.  The  computed  concentration 
g  r  c  tiles  ir.  the  x  direction  at  different  y 
locuticr.s  show  a  well -mixed  core  region  surrounded 
b\  concent  rat  ion  variations  in  the  outer  regions 
and  the  arms  of  the  vortex  for  Schmidt  number  of 
unity  The  character  of  the  concentration 


profiles  for  large  Schmidt  numbers,  typical  of 
liquids,  is  considerably  different  showing  sharp 
variations  of  concentration  between  adjacent 

stratifications  within  the  vortex  which  have  not 
diffused  because  of  the  low  molecular 

diffusivity  The  pdfs  constructed  from  such 

concent  rat i on  profiles  are  shouTi  in  Figures  2.  3, 

and  u  to  delineate  the  effects  of  vortex  strength, 
S c  r.m :  d  t  numb e  r  and  time. 

Figure  2  shows  the  pdf  as  a  function  of  the 
concentration  and  normalized  transverse  distance 
for  two  vortex  strengths  at  a  given  time 

and  Sc  -  u/V  of  uni  tv.  The  length  scale  chosen  to 
normalize  the  transverse  distance  approximately 
accounts  for  the  growth  of  the  vortex  with  time 
and  its  dependence  on  the  vortex  strength. 

Therefore ,  the  normalized  coordinates  X  - 

x/jUvt  Re  and  Y  -  y/JUi/X.  Re  are  convenient  in 
presentation  of  results.  The  probability  density 
function  shows  two  peaks  at  the  unmixed 
concentrations  outside  the  mixing  region. 
Probabilities  at  intermediate  concentrat ions 
develop  as  a  result  of  mixing  by  molecular 

diffusion  and  vortex  interaction.  As  the  vortex 
strength  increases,  the  probability  of  the 
concentrations  in  the  mid-range  (0.4  -  0.6) 

increases  in  the  core  of  the  vortex.  The  peaks 
at  Y  >  0  for  C  >  0.5  and  those  at  Y  <  0  for  C  < 
0.5  are  results  of  the  vortex  motion  bringing  the 
high  concentrations  from  lower  plane  to  the  upper 
and  vice  versa  The  counterclockwise  vortex 

motion  is  also  apparent  in  these  pdfs.  The 
continuing  line  of  peaks  on  each  side  joining  to 
the  peaks  at  C  -  0.0  and  C  -  1.0  result  from 
diffusive  regions  around  the  interfaces.  The 
shapes  of  the  pdfs  are  similar  to  those  computed 

in  shear  layers.  It  also  resembles  the 

qualitative  features  of  the  experimentally 

determined  pdfs  in  shear  layers. 

Figure  3a  and  2b  shows  the  time  evolution  of 
pdf  at  a  vortex  strength  of  Re  -  50.  At  early 
time  instant,  the  probabilities  of  intermediate 
concent i  at i ons  are  small  and  in  fact  the  shapes  of 
the  probability  peaks  at  the  minimum  and  maximum 
concent  rat  ions  are  similar  to  those  observed  for 
the  high  Schmidt  number  case  shown  In  Figure  4. 
At  a  later  instant.  the  effect  of  the  vortex 
becomes  more  pronounced  and  the  higher 
probabilities  of  mixing  appear  in  the  center  of 
the  vortex . 

Figure  4  shows  the  probability  densities  at 
large  Schmidt  numbers  typical  of  liquids  These 
pdfs  are  characterized  by  two  peaks  at  the  unmixed 
concentrations  with  a  minute  amount  of  mixing  in 
the  core  region.  The  shapes  of  the  probability 
peaks  can  be  approximately  correlated  by  the 
length  scale  chosen  in  the  transverse  direction 
These  pdfs  also  resemble  the  earlier  time  pdf  for 
Schmidt  number  of  unity 

One  convenient  measure  of  molecular  mixing 
can  be  defined  in  terms  of  the  probabilities  and 
concentrations  as 

f  ■  ?•  -Ti"  11  ■  4  So  P<C-V>(C  ■  0  ‘>)2dC'dv  (1) 

In  this  expressior. .  P'C)  is  the  normalized 
probability  This  parameter  car.  bf  referred  to  as 
"m.  i  xedne  s  s "  and  its  values  range  from  zero  fo4  no 


mixing  (C  -  0  or  1  only)  to  unity  for  complete 
mixing  (C  -  0.5’  everywhere).  The  values  of  this 
parameter  were  computed  for  a  range  of  values  of 
vortex  strengths  (Re  -  10-500)  and  Schmidt  number 
(Sc  -  1  -1000).  In  these  calculations,  2L  was 

taken  to  cover  the  whole  extent  of  the  largest 
vortex  considered  here  (L  -  10.0).  Empirical 

correlation  of  the  computed  mixing  factors 

suggests  an  expression  given  by, 

1  /2 

f  -  0.8  Sc'1/2  [^]  +  0.1  Re  (2) 

L  L 

The  "mixedness"  increases  continuously  with 

increasing  vortex  strength.  The  mixing  at  zero 
Reynolds  number  limit,  the  first  term,  corresponds 
to  a  diffusion  layer  growing  in  the  y  direction 
with  time.  The  second  term  is  the  mixing 
enhancement  by  the  vortex.  Current  results 

indicate  that  the  enhancement  factor  is  linear  in 
Re.  However,  there  appears  to  be  negligible 
dependence  of  this  factor  on  Schmidt  number. 

Figures  5a  and  5b  show  the  concentration 
distributions  in  the  (X,  Y)  plane  for  nonreacting 
and  reacting  vortices  of  a  nondioensional  strength 
of  50.  The  nonreacting  vortex  is  charac ter ized  by 
a  well-mixed  core  around  C  —  0.5  and  the  diffusive 
regions  surrounding  the  core  and  the  vortex 
arms.  The  vortex  with  finite  rate  chemical 
reactions  has  a  considerably  different 
character.  As  expected,  the  reactive  species  has 
been  depleted  in  the  core  region.  Around  the 
core,  the  spatial  mass  fraction  gradients  appear 
to  be  steeper  than  the  nonreacting  vortex.  Below 
the  vortex  center,  a  tongue  of  unreacted  species 
is  observed.  This  region  corresponds  to  the  zone 
between  the  outer  arm  of  the  vortex  and  the  core 
where  one  reactant  is  abundant. 

The  corresponding  pdfs  are  shown  in  Figures 
3a  and  3b.  Figure  3a  displays  the  typical 
features  of  nonreacting  vortex  pdf  at  early 
times.  In  comparing  the  two  pdfs,  the  reacting 
vortex  pdf  does  not  contain  the  vortex  core  peaks 
and  the  effect  of  the  reacted  core  on  the  pdf 
appears  as  modification  of  the  probability  peak 
topography  at  C  -  0.0.  The  new  probability  peak 
at  C  -  0.0  near  the  vortex  center  corresponds  to 
the  reacted  core  region. 
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